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Abstract. - A first-principle continuum-mechanics expression for the rolling friction coefficient 
is obtained for the rolling motion of a viscoelastic sphere on a hard plane. It relates the friction 
coefficient to the viscous and elastic constants of the sphere material. The relation obtained 
refers to the case when the deformation of the sphere ^ is small, the velocity of the sphere V is 
much less than the speed of sound in the material and when the characteristic time £^/V is much 
larger than the dissipative relaxation times of the viscoelastic material. To our knowledge this 
is the first "first-principle" expression of the rolling friction coefficient which does not contain 
empirical parameters. 



Rolling friction is one of the basic phenomenon men encounters in his everyday life since 
antediluvian times when the wheel was invented. The phenomenon of rolling friction has been 
interesting to scientist for a long time. Scientific publications on this subject range back to (at 
least) 1785 when Vince described systematic experiments to determine the nature of friction 
laws [0, and important scientists dealt with this problem among them O. Reynolds (|. The 
rolling friction is of great importance in engineering and science. From the speed of landslides 
Huang and Wang |Q argued that the rolling friction plays an important role even in geological 
processes, for a theoretical consideration see [Q. For its major importance the phenomenon 
has been studied intensively by engineers and physicists (e.g. [||), however, surprisingly few 
is known about its basic mechanisms. To our knowledge there is still no "first-principle" 
expression for the rolling friction coefficient available which relates this coefficient only to the 
material constants of the rolling body and does not contain empirical parameters. 

It has been shown that surface effects like adhesion |^ , electrostatic interaction Q , and other 
surface properties might influence the value of the rolling friction coefficient. Theoretically 
this problem was studied in Ref. ||^ where the authors propose a model of a surface with 
asperities to mimic friction (see also ||l^). In other studies [|ll|, |l^ it was argued that for 
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viscoelastic materials "the rolling friction is due very little to surface interactions: the major 
part is due to deformation losses within the bulk of the material" [jll] . Based on this concept 
the rolling friction coefficient was calculated in |Q where the deformation in the bulk was 
assumed to be completely plastic; then an empirical coefficient was introduced to account for 
the retarded recover of the material. 

In the present letter we also consider the rolling friction as a phenomenon appearing due 
to viscous processes in the bulk. We assume that energy losses due to surface effects may be 
neglected, compared to the viscous dissipation in the bulk. Thus we attribute the effect of 
rolling friction to viscous dissipation in the material due to time-dependent deformation. We 
assume that the only role of surface forces is to keep the rolling body from sliding. We use a 
quasi-static approach and obtain the rolling friction coefficient. 




Fig. 1. - Sketch of the rolling sphere, the right figure shows the appearance of the friction momentum: 
in the front the elastic stress is enhanced by the dissipative stress, while in the rear it is diminished. 



Fig. shows a sketch of the sphere of radius R and mass m rolling with angular velocity 
Vl which is related to its linear velocity VasV — Q.R. The center of the coordinate system 
OXY Z coincides with the center of the sphere (moving coordinate system), so that the plane 
is normal to the Z-axis and located at Zp < 0. Within our model we assume that the surface 
is much harder than the sphere, hence, the surface remains flat in the contact area while the 
sphere is deformed. Thus the contact area is a circle of radius a. The compression ^ = i? — jzpl 
of the sphere causes the displacement field in the bulk of the material u(f, t) where r gives 
the location of the point within the sphere at time t. The displacement field gives rise to a 
(tensorial) field of deformation 'V o u, and thus to the elastic part of the stress tensor aei ffj] 



Vou + Mov|-i/V-M 



EoIV ■ u. 



(1) 



Here, / is the unit tensor and Ei — and E2 — 



denote the elastic material 



3(l-2i/) 

constants with Y and v being the Young modulus and the Poisson ratio, respectively. Since 
the deformation depends on time via the time-dependent field u(r , t) , the dissipative part of 
the stress tensor reads 



V o u + u o V 



/V • u 



(2) 
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where u denotes the time derivative of the displacement field and 771 and 772 are the viscous 
constants of the sphere material ]T^ . The total stress 17 is a sum of both parts, a = a^i + o'dis- 
The force normal to the plane which acts on the sphere is caused by its own weight mg. It 
reads 

F = Ja-nds^nJJ dxdy {a^f + afj = riF^ , (3) 

where n is the (positively directed) unit normal to the plane, a^f and cr^^^ are the corresponding 
components of the stress tensors. The integration is to be performed over the contact circle 
Sa on the plane z — Zp. Correspondingly, the torque acting on the sphere with respect to the 
point O' at rest (no slipping conditions) is given by 

M = / {f - fo') x ■ n) ds , (4) 

where fo' = (Oj Oj ^p)- The torque M thus characterizes the rolling friction. 

Let the sphere move in y-direction, V — (0, V, 0), then the angular velocity and the friction 
torque are directed along the x-axis, i.e. d = (—£7,0,0) and M = (M, 0,0) (see Fig. For 
this geometry from Eq. (^) follows 

A/ = y J dxdy ya^^ {x,y,z ^ Zp) (5) 

The integration in Eq. (^ again is performed over the contact circle. 

To perform the calculation of the friction torque one needs the displacement field u^f, t) 
and the field of the displacement velocity u{f, t) . Generally, to find these quantities is a rather 
complicated problem. Nevertheless, one can use a quasi-static approximation, provided that 
the displacement velocities in the bulk are much smaller than the speed of sound in the sphere 
material and provided the characteristic time of the process of interest, estimated by r = £,/V , 



is much larger than the dissipative relaxation times of the material (see |13 for details). This 
allows to employ for the displacement field u the results of the static problem, i.e. of the 
Hertz contact problem which refers to a "slow" elastic collision of two spheres . In the case 
considered here one of the spheres (the plane) should have infinite radius. Let Ueiiv) be the 
solution of this (static) contact problem. It corresponds to the displacement field in a sphere 
under the compression ^ (see Fig. 0) if it would be at rest. In the quasi-static approximation 
the distribution Ueiiv) persists (i.e. it does not change with time) in the coordinate system 
OXYZ^ moving with the velocity V (with no rotation). In a body- fixed coordinate system 
OXYZ (see Fig. |^), which rotates with the angular velocity Vl and coincides at some time 
instant with OXYZ the displacements distribution is the same as in OXYZ ( i.e. Uei{r) ), 
while the displacement velocity distribution follows from the kinematic equation 

u(f) = (d • f X V)we; (r) = - ft {y dz - z dy) u^i (r) (6) 

which relates both coordinate systems; here 9x — dy — and dz — -§^- The elastic part 
of the stress tensor is then obtained by substituting Uei {r) into Eq. (|l|) . This part of the stress 
tensor corresponds to the static case when no torque acts on the sphere. Therefore, in the 
quasi-static approximation only the dissipative part of the stress tensor a j| should be used in 
Eq. (^) for the friction torque. Using kinematic Eq. ^ and definitions (^ and (||) one finds 
for the dissipative part 

''dts = -^{iydz-zdy)a^J{E,^r^,;E2^m) + 

+2al^{E, ^ irj2 - r/i/3);i?2 ^ m) " i^m + r/i/3)a,<J . (7) 
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where vl^i = ("ei' '"ep "ei) ^1 ^ ^2 ^ '72i etc. in the arguments of the elastic stress 
tensor denote that the elastic constants Ei, E2, etc. in Eq. (|^) should be replaced by the 
dissipative constants rji, 772, etc. 

For the case of small compression ^/i? <C 1 which is addressed here we estimate the 
magnitude of different terms in Eq. (|^) taken at the plane z = Zp fa —R, since these values 
of the stress tensor determine the friction torque (see Eq.(^). For the diagonal component of 
the stress tensor one has the known solution for the Hertz contact problem 

a'Jix, y, zp) = E^d^uli + {E2 - i?i/3) (a.<; + + d^^uli) 

3 I 

where a = b = \/2^R is the radius of the contact circle and F^i is the total elastic force, acting 
by the surface (in normal direction) on the sphere: 

F^AO^lj^f'H"' . (9) 

From Eq. (||) follows aff ~ at z = Zp. Applying geometric consideration of the 

deformation of the sphere one gets u^; ~ u^^ ^ ^^^'^/R^^^, and u^; ~ ^. The estimates 
dx dy ^ 1/a ^ {^R)~^/^ and dz ~ {S,R)~^^'^ may be obtained too. The expression for 
dz follows from the estimates for u^i, crff and the definition of crff. Finally we find for the 
deformation 

d.uli ~ dyu:^ ^ dzuli ^ ^/R « i^/Rf' 
d.ul ~ dyul ~ dzul ^ i/R « i^/Rf' 

dxul,^dyuli^dzuli^{^/Rf\ (10) 

Eqs. ( p^ ) further yield a^J ^ (f/i?)^^^. Using these estimates one can calculate the relative 
magnitude of the terms in Eq. (|^) and then compare their contribution to the rolling friction 
torque M. After these calculations one arrives at 



M = -nR j j dxdy ydy a^f {Ei <^ rii. E2 ^ m) 



(11) 



1/2 

where all other terms which are by a factor {£,/R) smaller are omitted. Integrating by parts 
the right-hand side of Eq. (|ll|) and taking into account that the stress tensor vanishes on the 
boundary of the contact circle one finally finds 



J dxdy alf {El ^m;E: 



M^VtRl / / dxdy alf {Ei ^ E2 <-> m) > (12) 



The expression in the curly brackets in Eq. ([l^ ) would coincide with the normal elastic force 
(see Eq. (||)) if it would contain the elastic constants Ei and E2 instead of the dissipative 
constants 771 and rj2. To perform the integration in Eq. (|l^ ) we apply the rescaling procedure 
proposed in [hsl. The coordinates are transformed due to a; = ax' , y — ay' and z = z' with 



m - lm\ fE2 + lEi 



(13) 



JI2 + Im) \E2 - \Ei 

so that the radius of the contact area in new coordinates is a — 6 = aa' = ah' . One obtains 



rlf{E,^^,-E2^m) - /^^i^+l^^^-Df^" 



dwi, dul 



\ dx dy dz 
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dz' \ 3 / \ dx' dy' dz' 



2 TTa'h' 



,/2 



where /3 = ^ ^^~|^\ As it follows from Eq. (^ after the appropriate rescaling the desired 
expression for the stress tensor with the "replaced" constants may be written in the same way 
as in (|^). Integration in Eq. ( p^ ) yields the final result for the torque of the rolling friction in 
the quasi-static regime 

M = k^^flRF^ . (15) 

It may be shown that in the quasi-static regime the normal force which acts on the sphere 
coincides with F^ of the static problem. Fig. |l| explains the appearance the friction torque in 
the rolling motion. The rolling friction coefficient finally reads 



^rol 



1 (3??2 -Vi f 
3 (3ry2 + 2?7i) 



(1-Z.2) (i_2^) 



(16) 



The latter expression relates the rolling friction coefficient to the elastic and viscous constants 
of the rolling body material. 

It is interesting to note that the rolling friction coefficient k^^i coincides with the material 
constant A in the dissipative force F^^^ — pA^S^, p = jj^-^VR for two colliding spheres 
(for derivation see pq]). Using this expression one can calculate the normal coefficient of 



restitution e for two spheres which collide with relative velocity f |13 . This coefficient which 
measures the energy loss upon collision reads 1 — e — Ci Ap^^^v^l'^ + • • • (C\ = 1.15344) [ p^ . 
Therefore, one can find the rolling friction coefficient from the restitution coefficient of the 
same material. To get an estimate of k^^i we take 1 — e — 0.1 for steel particles of mass 
m = 10^2 kg, colliding with impact velocity v = Ira/ s to obtain p^/^ — 10^ {kg^ /ms^y^^ and 
^rol — 10~^ s. Up to our knowledge for the case of a soft sphere rolling on a hard plane the 
dependence of the rolling friction coefficient on the rolling body velocity has not been studied 
experimentally. 

In conclusion, we found for the first time, to our knowledge, the first-principle expression 
which relates the rolling friction coefficient to the elastic and viscous properties of the material. 
Our calculations refer to the rolling motion of a viscoelastic sphere on a hard plane without 
slipping in the regime when the velocity of the sphere is much less then the speed of sound 
in the material and when the characteristic time of the process S,/V is much larger than the 
dissipative relaxation times of the viscoelastic material. 
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by Deutsche Forschungsgemeinschaft (Po 472/3-2). 
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